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Direct tissue imaging mass spectrometry (IMS) by matrix-assisted laser desorption 
ionization and time-of-f light (MALDI-TOF) mass spectrometry has become increasingly 
important in biology and medicine, because this technology can detect the relative 
abundance and spatial distribution of interesting proteins in tissues. Five thyroid cancer 
samples, along with normal tissue, were sliced and transferred onto conductive glass slides. 
After laser scanning by MALDI-TOF equipped with a smart beam laser, images were created 
for individual masses and proteins were classified at 200-|im spatial resolution. Based on 
the spatial distribution, region-specific proteins on a tumor lesion could be identified by 
protein extraction from tumor tissue and analysis using liquid chromatography with 
tandem mass spectrometry (LC-MS/MS). Using all the spectral data at each spot, various 
intensities of a specific peak were detected in the tumor and normal regions of the thyroid. 
Differences in the molecular weights of expressed proteins between tumor and normal 
regions were analyzed using unsupervised and supervised clustering. To verify the presence 
of discovered proteins through IMS, we identified ribosomal protein P2, which is specific 
for cancer. We have demonstrated the feasibility of IMS as a useful tool for the analysis of 
tissue sections, and identified the tumor-specific protein ribosomal protein P2. 

Keywords: Pathology; Proteins; Spectrometry, Mass, Matrix-Assisted Laser Desorption- 
lonization; Thyroid Gland; Neoplasms 



INTRODUCTION 

Papillary thyroid carcinomas (PTC) are the most common form 
of thyroid cancer, accounting for nearly 85% of primary thyroid 
malignancies. Commonly, the solitary nodule is a palpably dis- 
crete swelling within an otherwise apparently normal thyroid 
gland, and is incidentally detected by radiological evaluation 
(1). From a clinical standpoint, the possibility of neoplastic dis- 
ease is a major cause for concern in patients who present with 
thyroid nodules, because most conventional papillary carcino- 
mas present as asymptomatic thyroid nodules. Recent studies 
have demonstrated 2 molecular mechanisms that function in 
the carcinogenesis of thyroid cancer: the mitogen-activated 
protein kinase (MAPK) and phosphatidylinositol-3-kinase 
(PI3K)/Akt pathways (2, 3). Additionally, such pathways are re- 
lated to genetic alterations such as rearrangements in the "rear- 
ranged during transfection" {RET) or neurotrophic tyrosine ki- 



nase receptor 1 (NTRKl) genes, and the activation of point mu- 
tations in the BRAF gene (4). In a previous study, protein ex- 
pression patterns according to genetic variables have been 
demonstrated to correlate with the microscopic features, clini- 
cal manifestations, and prognostic characteristics of PTC (5). 
Although various protein alterations have been detected as di- 
agnostic markers or prognostic predictors, the challenge re- 
mains to identify the expression of new molecules in tumors. 

In biological research, protein expression profiling technolo- 
gy is a useful method to identify differential protein expression 
patterns and modifications. Analytical techniques with high 
sensitivity and increased throughput are required for the dis- 
covery of new biomarkers and new drugs. Recent progress in 
imaging mass spectrometry (IMS) has made it possible to iden- 
tify several cellular components such as proteins, drugs, and 
other endogenous molecules directly on tissue sections (6-10). 
IMS uses matrix-assisted laser desorption ionization time-of- 
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flight (MALDI-TOF) mass spectrometry (MS) technology. The 
matrix is applied on cryosectioned tissue and a very small area 
of the matrix- applied tissue is analyzed using MALDI-MS. Dif- 
ferences between the analyzed areas are displayed by the imag- 
ing program. Each analyzed area of the tissue section, which 
can be as small as > 200 jam, on a conductive surface such as 
gold-plated or indium tin oxide (ITO)-coated slide, is analyzed 
spot by spot using MALDI-MS after the application of the ma- 
trix. Using all the spectral data from each spot, the magnitudes 
of specific peaks on each spot can be displayed in terms of color 
intensity. In this way, spatial information on the tissue can be 
obtained. Based on the spatial distribution, region-specific pro- 
teins on a tissue can be identified by extracting proteins from 
the tissue, digesting them with trypsin, and analyzing the frag- 
ments using liquid chromatography with tandem MS (LC-MS/ 
MS). This technology using IMS has broad applications in the 
detection of new proteins in various fields. Compared with 
conventional imaging methods, the advantages of IMS are that 
it does not require the additional use of a specific antibody 
against the protein, and that it integrates histopathology and 
protein expression (11). For this reason, IMS provides superior 
information regarding distinct molecular arrangements in tis- 
sue sections. 

In the present study, we attempted to analyze thyroid sam- 
ples, including those from PTC, using IMS. Before analysis us- 
ing IMS, the diagnosis was confirmed by microscopic evalua- 
tion, immunohistochemical (IHC) staining of markers such as 
CK19, galectin-3, and RET, and pyrosequencing of the BRAF 
mutation, and then 5 homogenous samples of PTC were se- 



lected. After comparison of the differential molecular weight 
distribution between normal and tumor tissues using IMS, a 
protein of a specific molecular weight was identified by se- 
quencing combined with LC-MS/MS, and then the presence of 
this protein was reconfirmed by western blot analysis. 

MATERIALS AND METHODS 

Patient selection 

This study investigated 5 patients who had been diagnosed with 
papillary carcinoma in the thyroid and had undergone total 
thyroidectomy (Fig. 1). Four of the patients were female and 
one was male. Patient ages ranged from 48-68 yr (mean age, 
54.6 yr). Perithyroidal invasions of the tumor were observed in 
all patients, lymph node metastasis in 3 patients, and associa- 
tion with nodular hyperplasia in 2 patients. Tumor sizes ranged 
from 0.6-1.5 cm (average size, 0.96 cm). To prepare tumor sam- 
ples with identical molecular characteristics, the same configu- 
ration (BRAF mutation-positive, CK19 and galectin-3, and RET- 
negative), as confirmed by pyrosequencing and IHC staining, 
was selected. 

Tissue sections 

All samples were immersed in liquid nitrogen immediately af- 
ter excision. The frozen tissue was then sectioned into 10-|im 
thick slices at -20°C using a cryo-microtome (CM- 1850; Leica 
Microsystems Nussloch GmbH, Nussloch, Germany). The opti- 
mum cutting temperature (OCT) polymer can cause deteriora- 
tion of IMS signals; therefore, special care should be taken to 




Fig. 1. Gross and microscopic findings of tlie tliyroid. (A) Tlie riglit tliyroid sliowed a papillary carcinoma (black arrow) and nodular hyperplasia (black dash line arrow). Immuno- 
histochemistry was positive in CK19 (B) and Galectin-3 (C). 
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avoid OCT contamination in sections where it has been used to 
stabilize the tissue. The optimal procedure involves the use of 
OCT to adhere the tissue to the sample disc, without allowing it 
to come into contact with the sliced tissue. The frozen tissue 
sections were then thaw-mounted on ITO-coated glass slides 
(Bruker Daltonics Inc. Bremen, Germany). The prepared slides 
were dried in desiccators for about 20 min and then stored at 
-80°C until use. 

Sample preparation for imaging mass spectrometry 

Sinapinic acid (SA; Bruker Daltonics Inc.) was used as the pro- 
tein MALDI matrix and prepared as a 10 mg/mL solution in a 
50:50 solution of acetonitrile and 0.5% trifluoroacetic acid 
(TFA). An ImagePrep instrument (Bruker Daltonics Inc.) was 
used to spray a total of 3 mL of matrix solution on each tissue 
section. The optimal parameters (dry time, incubation time, 
and thickness) of the ImagePrep instrument were set to obtain 
a homogeneous matrix crystal on the tissue sections. After ma- 
trix application, the homogeneity of the matrix on the tissue 
was checked using the imaging function of Chip- 1000 (Shi- 
madzu Biotech, Kyoto, Japan). The ITO slide containing the tis- 
sue section was mounted on an MTP slide adapter (Bruker 
Daltonics Inc.), which was directly transferred to the MALDI 
mass spectrometer. 

Matrix-assisted laser desorption ionization and time-of- 
fliglit mass spectrometry 

IMS spectra were acquired in both positive- and negative-ion 
reflector modes with the aid of an Autoflex III MALDI-TOF 
mass spectrometer equipped with smart beam laser technolo- 
gy (Nd:YAG, 355 nm; Bruker Daltonics Inc.). MS data were ac- 
quired by averaging signals from 500 consecutive laser shots 
with a frequency of 200 Hz. All experiments conducted in this 
study were independently repeated 3 times. The spatial resolu- 
tion of each of the mass images shown is 200 |im. Prior to each 
data acquisition, external calibration was conducted using a 
mixed peptide calibration standard with a m/z range of 800- 
3,200. All obtained spectra were baseline-subtracted. 

Image analysis 

Data analysis software such as Fleximaging ver. 2.0, FlexAnaly- 
sis ver. 3.0 and ClinProTools ver. 2.1 from Bruker Daltonics (Bil- 
lerica, MA, USA) was employed to obtain the density maps of 
the species of interest and to check differences in the intensity 
of each spectrum. We manually analyzed normal and cancer 
regions according to supervised classification. Additionally, 
unsupervised classification was performed with Principal 
Component Analysis (PCA) using the ClinProTools software, 
and the results of this analysis were displayed using the Flexim- 
aging software. 



Tissue lysis and high performance liquid chromatography 
separation 

After tissue sectioning for IMS, the remaining tissue masses 
were separated into normal and cancer regions and subjected 
to further protein identification. To remove any remaining OCT 
and blood, the tissue sample was washed with PBS. The T-PER 
tissue protein extraction reagent (Pierce/Thermo Fisher Scien- 
tific Inc. Rockford, IL, USA), with added protease and phospha- 
tase inhibitors for preventing protein degradation, was used as 
the lysis buffer. The suspended tissue samples were Dounce- 
homogenized at 200 rpm for 2 min on ice followed by sonica- 
tion (5 sec x 2). The final tissue lysates were centriftiged at 10,000 
X g for 5 min at 4°C to remove intact cells and cell debris. 

For the separation of target proteins from tissue lysates, se- 
quential chromatography using C4 and C8 columns was per- 
formed. For the C4 column, the initial protein amount was 3 
mg, and this was fractionated into 60 fractions. The fractions 
that contained proteins with the target m/z values were pooled 
and separated further using a C8 column. The resulting frac- 
tions were verified by MALDI-MS, and then subjected to pro- 
tein identification using LC-MS/MS. 

In-gel digestion and mass spectrometry analysis 

Fractions of interest were digested by a trypsin in-gel digestion 
procedure following separation by sodium dodecyl sulfate 
polyacrylamide gel electrophoresis (SDS-PAGE). Bands of in- 
terest were manually cut out, and the gel pieces were destained 
and washed prior to in-gel digestion. In-gel digestion was car- 
ried out with 250 ng/ |iL sequencing grade modified trypsin 
(Promega, Madison, WI, USA) in 50 mM NH4HCO4 buffer (pH 
8.0) at 37°C overnight. The digested peptides were extracted us- 
ing 5% formic acid in acetonitrile and cleared by centriftigation 
at 1,000 X g for 5 min. The supernatant was dried in a SpeedVac, 
and the resulting tryptic peptides were analyzed by LC-MS/MS 
with a ftised-silica microcapillary C18 column (75 \xm internal 
diameter; 100 mm length). LC was performed with the follow- 
ing linear gradient conditions: 0 min 3% B, 15 min 8% B, 23 min 
16% B, 50 min 48% B, and 52 min 90% B (buffer A: 0.1% formic 
acid in H2O; buffer B: 0.1% formic acid in acetonitrile). The col- 
umn was eluted with 90% solvent B for 10 min at the end of 
each run, and then re-equilibrated with the initial solvent con- 
dition (3% solvent B) for the next run. The flow rate was 250 nL/ 
min. The separated peptides were subsequently analyzed by 
linear ion trap mass spectrometry (LIT-MS) on an LTQ spec- 
trometer (Thermo Fisher Scientific Inc., San Jose, CA, USA). 
The electrospray voltage was set at 2.0 kv; and the threshold for 
switching from MS to MS/MS was 500. The normalized colli- 
sion energy for MS/MS was 35% of the main radio frequency 
amplitude, and the duration of activation was 30 ms. All spectra 
were acquired in data- dependent mode. Each full MS scan was 
followed by nine MS/MS scans, corresponding to the 9 most 
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intense peaks in the full MS scan. The repeat count of each peak 
for dynamic exclusion was 1, and the repeat duration was 30 sec. 
The dynamic exclusion duration was set at 180 sec and the ex- 
clusion mass width was ±1.5 Da. The list size of dynamic exclu- 
sion was 50. 

Database analysis 

The acquired LC-ESI-MS/MS fragment spectra were analyzed 
in the BioWorksBrowser™ (version Rev. 3.3.1 SPl, Thermo 
Fisher Scientific Inc., CA, USA) with the SEQUEST search en- 
gines against the National Center for Biotechnology Informa- 
tion (http://www.ncbi.nlm.nih.gov/) non-redundant human 
database. 

The SEQUEST criteria for identifying peptides included the 
following: Xcorr values greater than 1.8, 2.2, and 3.5 for +1, +2, 
and +3 charge state peptides, respectively, a delta Cn value > 0.1, 
a fixed modification of carbarmidomethylation at cysteine (+58 
Da), and a variable modification of oxidation at methionine 
(+16 Da). 

Ethics statement 

This study protocol was approved by the institutional review 
board of the Konloilc University Hospital (IRB No. KUH 1210010). 
Informed consent was confirmed by the board. 

RESULTS 

Tissue profiling by matrix-assisted laser desorption/ 
ionization imaging mass spectrometry 

Hematoxylin and eosin (HE)-stained slides allowed the visual- 
ization of tumor lesions as dark pinkish areas and normal tissue 
as light pinkish areas (Fig. 2). We obtained spectra from tumor 
and normal areas of the thyroid tissue during the scan (Supple- 
mentary Fig. 1). These spectra clearly showed multiple peaks, 
representing patterns of protein expression in tumor and nor- 
mal tissue. For the evaluation of a particular protein expression 
pattern, these spectra were analyzed using unsupervised and 
supervised clustering analyses, and significant m/z values were 



acquired that allowed the classification of tissues as tumor or 
normal samples (Supplementary Table 1). 

Spectrometric images of molecules with different m/z values 
on frozen sections revealed various imaging patterns. Images 
reconstructed from supervised- classification analysis and PCA 
of the acquired datasets clearly differentiated tumor and nor- 
mal regions. The classification and PCA data revealed that the 
molecular distributions match the HE-stained images of excis- 
ed thyroid tissue samples. Of five cases, 4 tissue samples dem- 
onstrated that peaks corresponding to m/z values of 11,651- 
11,667 exhibited higher intensities in tumor lesions, whereas 
the peak with an m/zvahxe of 9,268-9,284 consistently showed 
higher intensities in normal thyroid. Color was used to repre- 
sent ions, and color saturation was a function of the relative in- 
tensity (see color reference bar). 

Protein identification of tumor tissues 

In the above sections, we have described the usefulness of IMS 
for protein profiling in tumor and normal thyroid tissue sam- 
ples. To identify differentially expressed proteins, we extracted 
proteins via lysis of surgically excised tissue specimens from 5 
patients. For direct MALDI-TOF MS analysis, proteins extracted 
from tumor and normal tissue samples were separated by chro- 
matography on C4 and C8 columns. The fractions contained a 
protein with an m/z value of 11,661, consistent with one of the 
significant m/z values obtained from the above supervised and 
unsupervised clustering analyses of MALDI-IMS data (Supple- 
mentary Fig. 2). 

For further evaluation, the pooled fractions were digested by 
a tryptic in-gel digestion procedure following SDS-PAGE. The 
resulting tryptic peptides from the pooled fractions were se- 
quenced using LC-MS/MS. LC-MS/MS results were analyzed 
using the BioWorksBrowser™ (version Rev. 3.3.1 SPl, Thermo 
Fisher Scientific Inc., CA, USA) for the identification of peptide 
sequences (Fig. 3). 

Peptide profiles analyzed by BioWorksBrowser™ software 
were matched to 64 amino acid sequences (55.7%) within 3 
unique peptides of ribosomal protein P2 comprised of 115 
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Fig. 3. Protein identification (include mass spectrum for protein in HPLC fractions). 
Identification of proteins of differential expression was carried out by digesting of 
separated proteins with trypsin and sequencing by LC-MS/MS. MS/MS spectra 
showing a peptide sequences. YVASYLLAALGGNSSPSAK of fraction 50 (pooled 
fraction 4): Xc value is 2.420, ACn value is 0.261 (A), NIEDVIAQGIGK of fraction 50 
(pooled fraction 4): Xc value is 4.347, ACn value is 0.394 (B), LASVPAGGAVAVSAA- 
PGSAAPAAGSAPAAAEEK of fraction 50 (pooled fraction 4): Xc value is 4.308 (C), 
ACn value is 0.608 and tolerances of m/z of the peptide were + 2 and + 1 for full 
mass and MS/MS fragmentation, respectively. 



^'^-^\ -100% 



m/z 11661: detecting mass 




Fig. 4. In imaging mass spectrometry, the distribution of detecting protein (1 1 ,661 m/z) displayed high concentration and corresponded with concentration of ribosomal protein 
p2 (1 1 ,657 m/z) in tumor lesion (A mass: 4 m/z). 



amino acids. The correspondence between the identified pro- 
tein and ribosomal protein P2 was compatible with high confi- 
dence levels. The 3 unique matching peptide sequences had 
m/z values of 1,870.11 (amino acids 3-21; sequence, R.YVASYL- 
LAALGGNSSPSAK.D), 1,257.43 (amino acids 50-61; sequence, 
K.NIEDVIAQGIGK.L), and 2,776.08 (amino acids 62-94; se- 
quence, K.LASVPAGGAVAVSAAPGSAAPAAGSAPAAAEEK.K) 



(Supplementary Fig. 3). 

In addition, the MALDI-IMS data were re-analyzed to inves- 
tigate whether the identified ribosomal protein P2, with an m/z 
value of 11,657, showed differential distributions in the tumor 
lesions of another thyroid sample. The tumor tissue showed 
high concentration of detecting protein with an m/z value 
11,661, corresponding to ribosomal protein P2 as identified by 
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LC-MS/MS (Fig. 4). Additionally, the molecular distribution 
matched the HE-stained image of the surgically excised thyroid 
tissue. Therefore, IMS showed that the localization of the pro- 
tein was consistent with the location of the tumor lesion. 

DISCUSSION 

One of the most commonly used strategies in proteomics in- 
volves the identification of new proteins that are differentially 
expressed between normal and tumor tissues. The aims of can- 
cer proteomics are to define a more precise molecular classifi- 
cation of tumors than is accomplished by conventional classifi- 
cation, and to discover meaningful biomarkers for diagnosis, 
development of therapeutic strategies, and prognosis. Recently, 
IMS has been highlighted as a promising tool for the visualiza- 
tion and identification of biomolecules in tissue sections (12). 
In addition to conventional methods such as genomics or pro- 
teomics, IMS has been considered to be of high importance for 
the discovery of new molecules. In the present study, we dem- 
onstrated the classification of tissue samples into tumor and 
normal tissues by the analysis of differential protein molecular 
weight distributions in human PTC samples. Among the differ- 
entially expressed peaks, ribosomal protein P2 was identified 
by LC-MS/MS as a protein that is differentially expressed in the 
tumor area. Consequently, these findings suggest that IMS may 
contribute to the identification of various proteins in future 
studies. 

Profiling and imaging by IMS has been previously applied to 
tissue samples from various diseases, including lung, brain, 
breast, ovary, stomach and thyroid cancer (13-18). Molecular 
protein signatures represent a unique data set with which to 
classify and correlate cHnically relevant information and out- 
comes with changing molecular events during the progression 
and treatment of disease (19). The approach taken by most 
published studies are comparative proteomic analysis whereby 
mass spectral features (m/z peaks) are correlated with a variety 
of clinical data, such as therapeutic regimen and overall out- 
come. Some studies have demonstrated the identification of 
specific molecular changes associated with progression of dis- 
ease. For instance, molecular signatures fi:om gHoma, ovarian 
cancer (20, 21) and non-small-cell lung cancer can be used to 
classify these lesions according to histological grade and to pre- 
dict the patients' prognosis (13, 14). However, such significant 
molecular data have clinical limitations in predicting diagnosis 
and survival, because of the heterogeneous distributions of dif- 
ferent proteins in spite of the fact that each tumor sample has 
the same origin. 

Although clinical limitations still exist, IMS maybe helpful in 
the identification of unknown biological variables in various 
tumors. Therefore, this technique will provide a common dis- 
ease-wide approach that can be applied to search specific dis- 



ease entities. Previous data suggest that IMS might be superior 
to immunohistochemical staining in identifying the site of ori- 
gin for various tumors (22). Furthermore, the particular protein 
profiles identified by IMS could be useful to improve the selec- 
tion of drug agents during the development of a therapeutic 
strategy. 

Our IMS analyses of PTC sections reveal that the distribu- 
tions of differentially expressed proteins show a separation be- 
tween normal and tumor tissues based on discrete molecular 
weight. In particular, recently developed analytical methods, 
such as supervised clustering using manual analysis and unsu- 
pervised clustering using PCA, were conducted to improve the 
separation of distributional prototypes and to identify specific 
proteins in thyroid tissue. Subsequently, the presence of ribo- 
somal protein P2 was identified by LC-MS/MS. A pubHshed 
study reports that ribosomal protein P2 is associated with tu- 
mor growth and could be a potential molecular target for anti- 
sense therapy of human malignancies, but its function in PTC 
is still unknoAAHi (23). 

This study has several limitations. First, the small number of 
thyroid samples results in a limited ability to demonstrate an 
association between the detected protein and PTC. Second, in 
matching the m/z value of expected protein to that of ribosomal 
protein P2, the m/z value might have been altered due to tech- 
nological limitations such as insufficient mass resolution in the 
higher m/z region, although our m/z value was within the al- 
lowable error limit. Such technological limitations can make it 
difficult to ascertain our conclusions. Therefore, progressive 
technological improvement is essential to improve accuracy. 

Our study successfully showed the feasibility of using IMS as 
a useful tool for the analysis of tissue sections. This method may 
provide a new direction in the study of cancer, by identifying tu- 
mor-specific markers, and thereby leading to elucidation of the 
mechanism of tumorigenesis and the identification of thera- 
peutic targets. Additionally, protein profiles obtained fi:om IMS 
may allow the accurate prediction of tumor behavior, diagnosis, 
and prognosis, and uncover the etiologies underlying idiopath- 
ic disease. 
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Supplementary Fig. 1 . Matrix-assisted laser desorption/ionization mass spectra taken a different locations within a papillary carcinoma. Over 1 00 different peaks could be de- 
tected with some of them having a distinct spatial distribution in the tissue, top, tumor area; bottom, normal area. 
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Supplementary Table 1. Comparisons between tumor and normal tissues according to supervised and unsupervised cluster analyses 
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Note: standard error of single mass analysis: + 4D. It is in bold that significant values are concluded through consensus between supervised and unsupervised clustering. SC, 
Supervised clustering using manual analysis; UC, Unsupervised clustering using principle component analysis; N, normal tissue; T, tumor. 
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Supplementary Fig. 3. Peptide profiles analyzed by BioWorksBrowserTM software were matched to 55.7% of amino acid sequences within 3 unique peptides (red circle) of ri- 
bosomal protein P2. 
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